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ABSTRACT 



,S^ . We present observations of the rotational ortho-water ground transition, the two lowest para-water transitions, and the ground tran- 

{ZLh' sition of ionised ortho-water in the archetypal starburst galaxy M82, performed with the HIFI instrument on the Herschel Space 

Observatory. These observations are the first detections of the para-H20(lH-Ooo) (1113GHz) and ortho-H20^(lii-Ooo) (1115GHz) 
lines in an extragalactic source. All three water lines show different spectral line profiles, underlining the need for high spectral resolu- 
tion in interpreting line formation processes. Using the line shape of the para-H^OC 1 1 1 -Ooo) and ortho-HoO^C 1 1 1 -Ooo) absorption profile 



o 

-t— • . _ . „ . . _ ^- _... _ . -. _... 

^ ' in conjunction with high spatial resolution CO observations, we show that the (ionised) water absorption arises from a ~ 2000 pc^ 

' region within the HIFI beam located about ~50 pc east of the dynamical centre of the galaxy. This region does not coincide with any of 

the known line emission peaks that have been identified in other molecular tracers, with the exception of HCO. Our data suggest that 
water and ionised water within this region have high (up to 75%) area-covering factors of the underlying continuum. This indicates 

^ ■ that water is not associated with small, dense cores within the ISM of M82 but arises from a more widespread diffuse gas component. 

J^^- ■ Key words. Line: formation - Galaxies: ISM - Galaxies: individual: M82 - Infrared: galaxies - Submillimeter: galaxies 
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1. Introduction nearby universe (e.g. van der Werf etal. I2010l l. As part of the 

HEXGAL guaranteed time key program (PI Giisten), we are 
High-resolution spectroscopy of far-infrared and sub-millimetre surveying the low-level water lines in different nuclear envi- 
water lines is an important tool for studying the physical ronments, performing velocity-resolved spectroscopy with the 
and chemical properties of the interstellar medium (ISM). Heterodyne Instrument for the Far Infrared (HIFI, de Graauw 
Absorption by terrestrial atmospheric water vapour has ground- g^ ^j 2010). In this letter we report on our first observations to- 
based studies of water in extragalactic systems limited to radio ^^^ds the central region of the archetypal starburst galaxy M82. 
maser transitions (such as the famous 22 GHz water line) or to ^e adopt a distance of 3.9 Mpc (Sakai & Madore fT999] l. 



2. Observations and data reduction 



X 

j^ ' a few systems with significant redshift (e.g. Combes & Wiklind 
■ - - ' [1997 Cernicharo etal. ,2006 Menten etal. 2008). Earlier satel- 
lite missions, such as ODIN and SWAS, did not have enough 
collecting area to detect the relatively faint ground transitions of 
water in external galaxies. ISO and, more recendy, Spitzer have Using the HIFI instrument onboard Herschel, we have observed 
provided the first systematic st udies o f water in the far i nfrared the ground transitions of ortho and para water, o-H20(lio-loi) 
(IR) regime (e.g. Fischer et al l 19991 Gonzalez- Alfonso I2004J. and p-H20(lii-Ooo), as well as the p-H20(2o2-lii) line, to- 
These missions, however, did not cover the frequencies of the wards the centre of M82 (RA=09''55'"52'.22 Dec=69°40'46".9 
water ground transitions and other low-level water lines. J2000). Observations were carried out in fast-chopping dual- 
Only with the launch of the HerscheQ satellite, with its large beamswitch mode using a wobbler throw of 3' for all obser- 
collecting area, have these transitions become accessible in the vations. The wobbler frequency was 0.8, 1.4, and 2.0 Hz for 

the 557 GHz, 988 GHz, and 1113 GHz observations, respec- 

' Herschel is an ESA space observatory with science instruments pro- tively. Nodding was peif ormed every ~ 40 seconds. Calibration 

vided by European-led Principal Investigator consortia and with partic- was achieved through hot/cold absorber measurements every 20 

ipation from NASA. minutes. The data were recorded using the wide-band acousto- 
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optical spectrometer, consisting of four units with a bandwidth 
of 1 GHz each, covering the 4 GHz IF for each polarization with 
spectral resolution of 1 MHz. 

Data were reduced using the HIPI^ and CLASqj software 
packages. Spectra were calibrated using HIPE and then exported 
to CLASS format with the shortest possible pre-integration (typ- 
ically ~ 40 seconds). For each scan we combined the four sub- 
bands in each polarization to create a 4 GHz spectrum. From this 
spectrum we computed the underlying continuum using the line- 
free channels; we then subtracted first-order baselines from indi- 
vidual sub-bands. The baseline-subtracted sub-bands were again 
combined and the continuum level added. This results in a noise- 
weighted 4 GHz spectrum for each scan. These spectra were in- 
spected for remaining baseline instabilities, and scans with dis- 
torted baselines were omitted. For each line we inspected the 
co-added result in both polarizations (H and V) separately. The 
continuum level was found to agree better than 10% for each 
frequency. The integrated line intensities agree within 25% with- 
out significant differences of the line profiles (although the latter 
comparison is limited by the signal-to-noise ratio, in particular 
for the 557 GHz observations). In the following we therefore use 
the noise-weighted average of both polarizations, which yields 
effective on-source integration times for the three lines are 750, 
2250, and 1300 seconds for the o-H20(lio-loi), p-H20(lii-Ooo), 
and p-H20(2o2-l ii) lines, respectively. 

Since HIFI's calibration is still preliminary, we used the the- 
oretical predictions based on the Herschel's expected surface 
accuracy and the geometrical aperture size to convert the an- 
tenna temperatures to flux density (Kramer 120061) . This yields 
462 JyK', 467 JyK-\ and 470 JyK' at 557 GHz, 988 GHz, 
and 1 113 GHz, respectively. The final spectra are shown at a ve- 
locity resolution of 20, 10, and 25 kms"' in Fig.[T](top). 

3. Results 

All three water lines have been detected with high significance, 
demonstrating that faint (few mK) broad lines can be observed 
with HIFI. The o-H20(lio-loi) line is detected in emission and 
shows a double-peaked line profile. Within noise uncertainty, 
no emission (or absorption) is detected on the systemic veloc- 
ity of M82 (vlsr = 225 kms"'). Both components are well fit 
by Gaussian profiles, with line parameters given in Table [T] We 
detect a continuum level of 12.8 Jy beam"' at 557 GHz. 

Our 1 1 13 GHz spectrum shows two absorption features, one 
centered at the systemic velocity of M82 corresponding to the 
p-H20(lii-Ooo) line, and a second stronger feature blue-shifted 
by ~ 1.8 GHz with respect to the systemic velocity. We identify 
the blue-shifted absorption feature as o-H20^(lii-Ooo) (Viest - 
1115.186GHz, Murtz etal. 1998). Both absorption line profiles 
are identical within the uncertainties (see Fig.|2]). The absorption 
profile is approximated reasonably well by a single Gaussian 
(see Table [T] for the line parameters). The continuum level de- 
tected at 1 1 13 GHz is 78 Jy beam"' . 

The p-H20(2o2-lii) line at 988 GHz is detected in emis- 
sion. Its line profile differs from both the o-H20(lio-loi) emis- 
sion profile and thep-H20(lii-Oo()) absorption profiles. It shows 
emission between Vlsr = 50 - 500 kms"' (similar to the ve- 
locity range covered by CO at the same spatial resolution, see 
below) and is almost flat-topped for velocities between 100 and 
400 kms"'. Thus the spectrum does not indicate absorption (or 
absence of emission) at the systemic velocity. The spectrum 



^ Herschel Interactive Processing Environment 
^ http://www.iram.fr/IRAMFR/GILDAS 



can be decomposed into two Gaussian profiles, with parameters 
given in Table[Tl The continuum flux detected at 988 GHz is 61 

Jybeam"'. 



4. Discussion 

4. 1 . The line profiles 

Given the large body of high spatial resolution observations of 
molecular gas tracers published for M82, the line profile of the 
water lines can be compared to other data to learn more about the 
location and extent of the water emitting/absorbing regions in the 
disk. We here compare the water line profiles to the high spatial 
resol ution ( 3.5") '^CO(7=1^0) data cube obtained by Walter 
etal. ( I2OO2I 1: '^CO (CO thereafter) is the best-studied tracer of 
the molecular gas in M82. We first compare the CO spectra in 
beams synthesized to the same spatial resolutions as the HIFI 
beams. From the comparison of the p-H20(lii-0()o) absorption 
spectrum to CO, it is apparent that the absorption is not only de- 
tected in the pronounced absorption feature close to the systemic 
velocity, but also at velocities in the wings of the CO profile (Fig. 
[T]bottom). It is therefore tempting to speculate that the lack of 
absorption at certain velocities has a geometrical origin, i.e., that 
gas at these velocities is located behind the continuum. This is 
also supported by the shape of the p-H20(2o2- In) line emission 
profile, which shows that water is abundant in the gas phase of 
M82 at all velocities where CO (i.e. molecular gas) is present. 
The very good correspondence of the o-H20^(lii-Ooo) absorp- 
tion profile near the systemic velocity suggests that the ionised 
water traces the same gas as is detected in the water absorption. 
A closer inspection of the o-H20^(lii-0()o) profile shows, how- 
ever, a lack of absorption in the red wing of the line profile (see 
Fig|2]i. The blue wing is only partly covered by our spectrum and 
shows emission at the very blue edge. This could come from cal- 
ibration uncertainties at the edge of the IF-band, but it is unclear 
whether this feature is an artifact or is real. 
Our finding that the p-H20(lii-Ooo) line is observed in absorp- 
tion while the p-H20(2o2-lii) line is detected in emission can 
be used to obtain an estimate of the excitation temperatures of 
both lines. We used the dust model by Siebenmorgen & Kriigel 
(I2007I I to estimate a background temperature of 1 8 K and 20 K at 
988 and 1 1 13 GHz, respectively. This implies Tex < 20 K for p- 
H20(lii-0oo) and Tex > 18Kforp-H20(2o2-lii) (or Tex ^ 19K 
if both lines are close to LTE). Given the complexity of the water 
energy level diagram and the various level population channels 
(collisional or radiative), detailed models will be required for in- 
vestigating the underlying excitation mechanisms. 

Owing to the much larger beam size of the o-H20(lio- 
loi) observations, a comparison to the other water lines is not 
straightforward. In comparison to CO, however, it is apparent 
that the o-H20(lio-loi) emission arises exclusively from veloc- 
ities in the line wings of the CO spectrum. These velocities 
mainly correspond to emission in the southwestern and north- 
eastern molecular lobes that are located within the 41" beam (see 
Fig|3]for the CO distribution compared to the HIFI beams), and 
yet we cannot rule out that the lack of o-H20(l lo-loi) emission 
near the systemic velocities is due to self absorption in the water 
line profile. Observations of the water lines within the molecu- 
lar lobes will be required to obtain similar spatial coverage of 
all lines, thus allowing detailed analysis of different line profiles 
and modeling of the water excitation. These observations have 
been approved and will be presented in a forthcoming paper 

We retrieved sub-beam location information by comparing 
CO line profiles at the highest spatial resolution (3.5") to the p- 
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Fig. 1. Top: Spectra of the o-H20(l lo-loi), p-H20(l n-Ooo) and p-H20(2o2-l 11) lines towards the centre of M82, with Gaussian-fit 
profiles superposed. The line parameters are given in Table [1] Bottom: Water-line profiles superposed on the €0(7=1^0) line 
profiles within beams synthesized with the same spatial resolution as the water data. The velocity scale is relative to the systemic 
velocity of M82 (vlsr = 225 kms"'). The intensities of the CO line profiles have been normalized to unity, and the water line 
profiles scaled to provide close matches to the CO profiles in the line wings. The CO profile in the middle panel has been inverted 
for better comparison to the water absorption profile. 

Table 1. Line parameters derived from a Gaussian fit. 



Line 



Viest beam size 

[GHz] ["] 



T^ peak S^ peak 

[mK] [Jybeam-'] 



L Vlsr dv S, cont. 

[Jy beam"' kms"'] [kms"'] [kms"'] [Jybeam"'] 



o-H20(lio-loi) 556.936 

p-H20(2o2-lii) 987.927 

p-H20(lii-Ooo) 1113.343 

o-H20+(lii-Ooo) 1115.186 



41 


4.5 ±1.1 


2.1 ±0.5 


200 ± 33 


110±8 


90 ±15 


12.8 (540 //m) 




5.0 ±1.1 


2.3 ±0.5 


170 ± 29 


308 ±6 


68 ± 11 




23 


8.9 ±2.0 


4.2 ±0.9 


273 ± 90 


100 ±9 


61 ±20 


60.0 (300 /im) 




9.3 ±2.0 


4.4 ±0.9 


1005 ± 130 


250 ± 15 


220 ± 27 




20 


-12.0 ±1.6 


-5.6 ±0.7 


-432 ± 28 


240 ±2 


72 ±6 


78.1 (270 //m) 


20 


-21.4 ±1.6 


-10.1 ±0.7 


-836 ± 29 


242 ±2 


77 ±5 





H20(lii-0oo) and o-H20^(lii-0()o) line profiles. The results are 
displayed in Fig. [3] where we show selected CO spectra overlaid 
on the p-H20(lii-Ooo) absorption profile. We find that the CO 
line profile is in very good agreement with the water absorption 
profile in only a small region within the HIFI beam. The region 
delineated in Fig. [3] corresponds to a small strip orthogonal to 
the molecular disk of M82. The CO profiles east and west of this 
region show significant shifts of their line centroids compared to 
the peak of the water absorption due to the rotation of the molec- 
ular disk. Interestingly, the CO line profiles that match the shape 
of the main water absorption profile also show a blue wing. This 
CO emission corresponds to molecular gas in the outflow of M82 
(Walter etal. 2002). This shows that the water absorption in the 
blue line wing arises not only from the northwestern molecular 



lobe (which has similar velocities and is at least partly covered 
by the HIFI beam), but also from gas in the outflow. 



4.2. The origin of tlie H2O/H2O+ absorption 

We used the observed absorption depth to investigate whether 
the p-H20(lii-0()o) and o-H20^(lii-Ooo) absorption arises from 
small cores or from a more widespread gas phase of the ISM. 
For this we assume that both lines are optically thick. This im- 
plies that the absorption depth directly measures the continuum 
covering factor, once the continuum distribution is known. We 
here use the 3 mm continuum shown in Fig|4] as a proxy for 
the submm continuum distribution. Although the continuum at 
this wavelength is not automatically a good indicator of the 



A. WeiB et al.: HIFI spectroscopy of low-level water transitions in M82 




200 



100 -100 

Velocity (km/s) 



-200 



Fig. 2. Line profiles of the p-H20(lii-Ooo) (black) and the o- 
H20^(lii-0oo) (red) lines. The velocity scale is relative to the 
systemic velocity of M82. Both profiles have been normalized 
to a maximum absorption of unity. 



thermal dust emission, because of possible free-free contami- 
nation, the spatial similarity between the the 3 mm and radio 
continua at high spatial resolution, as well as the similarity of 
the 3 mm if smoothed to similar resolution to the published 
350 /im maps (Leeuw & Robson .2009j . justifies this approach. 
The total continuum emission of M82 at 1113 GHz was esti- 
mated using the dust model by Siebenmorgen & Kriigel (120071) . 
Combining this with the observed 3 mm continuum distribution, 
we estimate S'lincHz = 13.5 Jy within the region of the water 
absorption. This translates into a continuum covering factor of 
40% forp-H20(lii-Ooo) and 75% foro-HzO+On-Ooo). This ap- 
proach yields predicted continuum fluxes within the HIFI beams 
that agree within 10% or better with the observed continuum 
fluxes listed in Table[T] at all three frequencies. Since the wa- 
ter absorption region covers ~ 2000 pc^, these numbers strongly 
suggest that the p-H20(lii-Ooo) ando-H20^(lii-Ooo) absorption 
arise from a widespread gas phase of the ISM and not from 
dense cores (which are expected to cover a much smaller frac- 
tion of the continuum).. Considering the high critical densities 
for coUisional excitation of water (a few times 10^ cm"^, Faure 
etal. 12007] ). this also implies that the water excitation cannot be 
driven by collisions, but is most likely dominated by the IR field. 
This conclusion agrees with findings in other IR bright galaxies 
such as Arp220 and Mrk231 (e.g. Gonzalez- Alfonso etal. 120041 
Gonzalez-Alfonso et al.'2010'). 

We can also use the absorption depth to calculated the ab- 
sorbing gas column densities assuming that the absorption is 
not intrinsically saturated and using the continuum strength de- 
rived above (13.5 Jy). Applying the method by Menten etal. 
(120081 we derive J Tappdv of 38.6 kms"' and 107.2 kms"' 
from the p-H20(lii-Ooo) and o-H20^(lii-Ooo) absorption pro- 
files. With the low-excitation temperature approximation, Ni - 
1^— J Tappdv, this yields ground-level column densities of 

A^(p-H2O000) = 9.0 ■ lO'^^cm-2 and A^(o-H20+ 000) = 2.2 ■ 
lO'^cm"^. For the latter, the fine structure splitting of o-H20^ 
(see e.g. Ossenkopf et al. '2010) has been taken into account. 
Assuming an ortho-to-para ratio of 3: 1 this results in lower limits 
for the total column densities of A^(H20) > 3.6 ■ 10'"* cm"^ and 
A^(H20^) > 2.9 ■ 10'"* cm"^. By converting the measured CO in- 




Fig. 3. HIFI beam sizes (FWHM) for our observing frequencies 
shown on the integrated CO(y=l— >0) distribution (Walter etal. 
2002^ in the central region of M82. The green cross indicates 
our pointing centre which corresponds to the dynamical centre 
of M 82. The red crosses indicate the locations where the CO 
profiles at 3.5" resolution match the p-H20(lii-Ooo) absorption 
profiles, and the blue crosses indicate locations where the CO 
profile significantly differs. The spectra to the top and bottom 
show example CO spectra superposed on the water absorption 
profiles. Their positions have also been marked by a circle for 
better visualization. 



tensities in the absorption region to proton column densities us- 
ing an Xco conversion factor as given in Walter et al. (120021) . we 
derive limits for the fractional abundances of [H20]> 4.0 • 10"' 
and [H20^]> 3.3-10"'. These numbers are 2-3 magnitudes lower 
than the water abundances derived for Arp220 and Mrk231 
(~ 10 '', Gonzalez-Alfonso etal. 2004, Gonzalez-Alfonso etal. 
2010), but approach values found in molecular cloud cores, the 
diffuse gas in the Milky Way, and in the z = 0.89 absorption 
system PKS 1830- 211 (e. g. SneU etal. l2000l van der Tak etal. 
120101 Menten et al. dSOQgl) . 

4.3. Comparison to other molecular gas tracers 

As shown in Fig.|3]the water absorption region is not associated 
with a particular CO-bright region in M82, but is located be- 
tween the northeastern molecular lobe and the central CO peak. 
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This indicates that the water absorption is not associated with 
the bulk of the molecular gas in M82 or that the underlying 
continuum distribution favors an absorption in regions that are 
not bright in CO. The latter can be ruled out by comparing the 
spatial distribution of the water absorption to the radio (Wills 
et al. 119971 ) and mm-continuum ( WeiB et al. 12001 al l as a proxy 
for the spatial distribution of the submm continuum (e.g. using 
the FIR-radio correlation). Neither tracer provides evidence that 
the submm continuum is in particular pronounced towards the 
region of the water absorption (see Fig.|4|, which suggests that 
the small size of the water absorption is related to the ISM prop- 
erties rather than to the continuum distribution. 

We have further compared the spatial distribution of the 
dense gas as traced by HCN (Brouillet & Schilke I1993I I and 
H'-^CO^ (Garcia-Burillo et al. i2002i) to the water absorption re- 
gion. These tracers, however, closely follow the distribution of 
the CO emission, which rules out any interpretation that the wa- 
ter absorption occurs predominately in high-density regions. 

A possible mechanism to efficiently release water from dust 
grains into the gas phase are shocks (e.g. Cernicharo et al. 11 9991 1 . 
A comparison to the SiO distribution in M82 (Garcia-Burillo 
et al. 1200 11 1 shows that there is no indication of strong shocks 
in this particular region. The SiO observations, however, do not 
rule out the possibility that weaker C-type shocks could release 
H2O efficiently from the dust grains into the gas phase without 
disrupting the grains themselves. In this context it is interesting 
to note that the water absorption is located near the expected lo- 
cation of the x2 orbits (and their intersection with the xl orbits) 
proposed in the context of gas motions in the stellar bar poten- 
tial of M 82 (e.g. Greve etal. 2002). Shocks associated with the 
orbital intersection would provide a plausible mechanism for re- 
leasing water into the gas phase, and the geometry of the orbits 
would naturally explain why the absorption is only detected to 
the east of the dynamical centre of M 82 (see Fig. 2 in Greve 
etal. l2002l i. 

On the other hand molecular abundances in M82 are thought 
to be largely influenced by chemistry in photon-dominated re- 
gions (PDRs, see e.g. Martin et al. i2006) . UV photons photo- 
dissociate water and PDR models predict that significant water 
gas phase abundances should only be present in well-shielded, 
high column density regions (see e.g. Meijerink & Spaans .2005) . 
This interpretation is challenged by our finding that the water 
absorption region in M82 is located in a relatively low column 
density region within the disk of M82. Furthermore the water 
absorption coincides with a peak of the spatial distribution of 
HCO (Garcia-Burillo etal. .2002. see Figgli. HCO is a tracer of 
PDRs, and it is the only molecule so far imaged at high spa- 
tial resolution that shows enhanced emission in this region. The 
similarity of the p-H20(lii-Ooo) ando-H20^(lii-0()o) absorption 
profiles shows that ionizing photons are present in the absorbing 
medium. Therefore UV dissociation could suppress the water 
gas phase abundances, leading to much lower water abundances 
in M82 compared to those derived in Arp 220 and Mrk 231. This 
is in line with studies of other molecules with low ionization 
potential, such as as NH3 and HNCO, for which photo dissocia- 
tion has been suggested as a main driver for the low abundances 
observed in M82 (WeiB etal. l2001bl Martin etal. l2009l l. 
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